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Abstract 


In  this  paper  we  evaluate  various  proposed  VLSI  models  of  computa¬ 
tion.  While  there  is  consensus  on  the  appraisal  of  chip  area,  controversy 
remains  with  regard  to  computation  time.  Thus  we  have  analyzed  in  detail 
the  propagation  of  signals  on  dispersive  lines.  The  results  are  expressed 
in  terms  of  adimensional  parameters  characteristic  of  any  given  fabrica¬ 
tion  technology.  The  conclusion  is  that  both  current  and  projected  silicon 
technologies  fall  within  the  realm  of  the  capacitive  model,  where  a  dis¬ 
persive  line  can  be  replaced  by  a  capacitance  proportional  to  its  length. 
Diffusion  phenomena  appear  therefore  to  exceed  the  present  VLSI  horizon. 
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A  CRITIQUE  AND  AN  APPRAISAL  OF  VLSI  MODELS  OF  COMPUTATION 


1 .  Introduction 

The  central  question  in  the  design  and  analysis  of  algorithms  is  the 
definition  of  the  model  of  computation  to  be  adopted.  Indeed,  “performance” 
becomes  meaningful  only  in  relation  to  a  given  model.  This  model  is 
normally  the  simplified  abstraction  of  a  class  of  real  or  imaginary  machines; 
for  example,  the  RAM  or  Random- Access-Machine ,  is  the  model  of  practically 
the  totality  of  existing  (Von  Neumann  type)  processors.  The  model  of 
computation  is  the  simplest  possible,  compatibly  with  the  requirement  of 
being  realistic.  In  other  words,  while  a  model  aims  at  capturing  the 
essential  traits  of  a  system  or  technology,  its  simplicity  is  what  enables 
theoretical  appraisals  of  performance.  . 

Very-Large-Scale- Integration  (VLSI),  as  a  computing  environment, is  no 
exception.  Indeed  considerable  attention  has  been  paid  [1][2][3][4] 
to  the  definition  of  a  suitable  model.  The  basic  parameters  of  any  VLSI 
computation  model  are  chip  area  A  and  computation  time  T.  VLSI  systems 
display  a  trade-off  between  these  two  parameters,  each  of  which  represents 
a  well-defined  cost  aspect:  chip  area  is  a  measure  of  fabrication  cost  and 
computation  time  is  a  measure  of  operating  cost. 

A  general  feature  of  all  proposed  —  and  presumably  of  all  future  —  VLSI 
models  of  computation  is  that  a  chip  is  viewed  as  a  computation  graph, 
whose  vertices  are  called  nodes  and  whose  arcs  are  called  wires .  Nodes  are, 
by  and  large,  devices  and  are  responsible  for  information  processing  (com¬ 
putations  of  boolean  functions);  wires  are  just  electrical  connections , and 
are  responsible  for  both  transfer  of  information  and  distribution  of  power 
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supply  and  timing  waveforms. 

A  given  computation  graph  is  to  be  laid-out  in  conformity  with  the 
rules  dictated  by  technology.  These  rules  are  geometric  constraints  on 
admissible  layouts  and  typically  concern  widths  of  wires  and  transistor 
regions,  clearances  between  wires,  transistors,  etc.,  number  of  metallic 
layers,  permissible  orientations,  etc..  Once  a  layout  —  that  is,  a  legal 
planar  embedding  of  the  computation  graph  —  has  been  produced,  the  chip 
area  A  is  normally  the  area  of  the  smallest  rectangle  inscribing  the  lay¬ 
out,  and  is  the  sum  of  the  areas  of  wires,  transistors,  and,  possibly,  of 
some  wasted  space.  More  formally  we  have: 

Area  Assumptions 

Al.  (Wire  area)  All  wires  have  minimum  width  X  >  0  (which  includes  both 

the  actual  wire  width  and  the  clearance  between  wire  and  any  other  chip 
region)  and  at  most  v  ^  2  wires  can  overlap  at  any  point  (hypothesis  of 
bounded  number  of  layers).  [All  models.] 

2 

A2.  (Transistor-port  area)  Transistors  and  I/O  ports  have  minimum  area  2:  \ 
[All  models.] 

2  2 

A2.1  Transistors  and  I/O  ports  have  fixed  area  c^X  and  CpX  , 

respectively,  for  constants  c  and  Cp  [Brent-Kung  [2];  Chazelle- 
Monier  [  4 ] ]  . 

A2.2  The  chip  is  subdivided  into  compact  regions,  called  "self-timed"; 
within  a  self-timed  region  A2 . 1  holds,  while  drivers  of  inter¬ 
region  wires  have  area  proportional  to  the  wire-length  [Thompson 
[  3 ] ;  Seitz  [  5 ] ] - 

A3.  (Chip  area)  The  chip  area  A  is  at  least  the  sum  of  the  area  of  the 

wires,  of  the  transistors,  and  of  the  I/O  ports,  and  it  is  at  most  the 
area  of  the  smallest  rectangle  (or  convex  region)  enclosing  a  legal 
layout  of  the  graph.  [All  models.] 
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These  rules  are  quite  simple  and  uncontroversial.  Indeed  no  difficulty  arises 
in  appraising  the  area  of  a  given  computation  graph. 

Radically  different  —  as  to  a  consensus  among  researchers  —  is  the 
situation  regarding  the  computation  time  T.  To  acquire  the  necessary 
perspective,  let  us  call  "an  elementary  action"  the  change  of  output  of  a 
transistor  and  the  transmission  of  this  change  on  the  wires  connected  to  this 
output.  Thus,  given  a  computation  graph  —  which  supports  a  prescribed 
algorithm  —  the  designer  can  describe  the  execution  of  the  algorithm  as  a 
sequence  of  sets  of  elementary  actions.  In  other  words,  execution  is  con¬ 
veniently  modeled  by  a  single-source/single-destination  (corresponding  to 
begin  and  end,  respectively)  directed  acyclic  graph,  whose  arcs  correspond 
to  elementary  actions.  Each  arc  is  weighted  with  the  time  taken  by  the 
action  it  represents.  This  knowledge,  in  principle,  seems  quite  adequate 
for  the  evaluation  of  T,  by  simply  taking  the  value  of  the  most  time-con¬ 
suming  source -destination  path  in  the  acyclic  graph.  The  difficulty  lies, 
however,  in  the  assignment  of  values  to  the  arc  weights.  Indeed,  the 
proposed  computational  models  basically  differ  in  this  weight  assignment. 

More  formally  we  have: 

Time  Assumptions 

Tl.  (Propagation  time  along  a  wire). 

Tl.l  A  bit  requires  a  constant  time  t  to  propagate  along  a  wire, 

irrespectively  of  its  length.  (Brent-Kung) .  (We  refer  to  this 
case  as  the  synchronous  model. ) 

Tl.2  A  bit  requires  a  time  O(logi)  to  propagate  along  a  wire  of  iength 
l  (Mead-Conway ;  Thompson).  (We  refer  to  this  case  as  the 
capacitive  model. ) 
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T1.3  A  bit  requires  a  time  0(4  )  to  propagate  along  a  wire  of  length 
l  (Seitz;  Chazelle-Monier) .  (We  refer  to  this  case  as  the 
diffusion  model. ) 

T2.  (Algorithm  time)  The  computation  time  of  an  algorithm  is  the  time  of 

the  longest  sequence  of  wire  propagation  times  between  beginning  and 

completion  of  the  computation.  [All  models.] 

The  choices  for  T1  reflect  the  profound  controversy  on  VLSI  computation 

time.  In  a  preliminary  analysis,  one  is  tempted  to  conclude  that  T1.2  is 

the  most  realistic  choice.  Indeed,  a  wire  is  characterized  by  a  resistance 

and  a  capacitance  which  (in  a  given  fabrication  technique)  both  grow 

linearly  with  the  wire  length;  therefore,  the  time  constant  of  the  transistor 

2 

load  grows  proportionally  to  l  ,  whence  the  conclusion  Tl.3.  Notice  that 
the  computational  implications  of  Tl.3  —  as  noted  by  Chazelle-Monier  in 
[4]  —  are  drastic.  Indeed,  chip  wires  of  substantially  different  lengths 
are  ruled  out  and  connections  must  exist  only  between  devices  in  very  close 
proximity.  As  a  consequence,  the  only  permissible  computation  graphs  are 
of  the  mesh  type  (or  closely  related),  which  rules  out  very  fast  parallel 
computation,  such  as  performed  by  computing  structures  of  the  type  of  the 
shuffle -exchange  [6],  the  cube-connected-cycles  [7],  or  the  tree-connected 
machine  [  8 ] . 

Asymptotically,  the  line  of  arguments  sketched  above  is  unimpeachable, 
and  therefore  —  for  the  theoretician  of  algorithmics  —  valid,  since 
asymptotic  analysis  is  the  cornerstone  of  concrete  computational  complexity. 
However,  the  asymptotics  of  VLSI  have  a  much  closer  horizon  than,  for 
example,  the  asymptotics  of  the  Turing  machine.  This  horizon,  in  fact,  is 
set  by  realistic  sounds  on  the  expectations  —  in  the  current  technology  — 
of  minimum  feature  size  and  maximum  chip  size. 
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Within  this  horizon,  the  line  parameters  must  be  weighed  against  the 
nonnegligible  output  impedance  of  the  driving  transistor  and  the  input 
impedance  of  the  driven  transistor.  To  appraise  this  interaction,  it  is 
therefore  appropriate  to  take  a  critical  look  at  the  actual  physical 
phenomena  occurring  during  an  "elementary  action". 


A  mathematical  model  of  wire  switching 


Perhaps  the  most  characteristic  feature  of  present-day  VLSI  technology 
is  the  fact  that,  irrespective  of  the  choice  of  the  devices  (MOS-FET  versus 
bipolar,  for  example)  wires  are  realized  as  dispersive  lines.  This  nature 
of  wires  is  what  determines  the  time  behavior  of  networks  (and  must  be 
reflected  in  the  computation  model)  and  the  choice  of  devices,  or  of  their 
operating  regimes,  has  a  nonessential  effect  on  it.  Therefore,  with 
reference  to  dispersive  line  VLSI  technology,  any  reasonable  device  selection 
is  representative  of  the  general  problem. 

In  particular,  we  shall  carry  out  our  analysis  with  reference  to  the 
CMOS  technology  [9].  In  figure  la  we  have  illustrated  the  circuit  being 
considered.  T^  is  an  n-channel  MOS  transistor,  initially  cut-off.  Its 
drain  load  —  that  is  the  wire  AB  and  the  gate  capacitance  of  the  driven 


transistor  T2  ~  is  initially  charged  to  voltage  Vq.  So,  with  reference  to 
(Ids»vds^  characteristic  curves  of  figure  lb,  P^  is  the  initial  operating 
point  of  T^.  At  t  =  0  a  step  voltage  v^  =  Vq  is  applied  at  the  gate  of  T^; 
after  a  time  Tq  —  negligible  with  respect  to  the  other  intervening  times  — 


the  current  Iq  corresponding  to  v^  =  Vq  is  established  and  the  operating 
point  moves  to  P,, .  From  this  point  on,  the  operating  point  moves  on  the 


v  =  Vq  curve  towards  the  origin  and  the  transistor  load  discharges  through 
8  — 

the  channel.  It  is  our  objective  to  analyze  this  phenomenon. 
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v  =  V, 


Figure  1. 


The  CMOS  configuration  and  the  transistor  characteristics. 

— i„t~ — — . - 


v(x,t) 


(a)  (b)  V 

Figure  2.  The  model  (a)  and  the  idealized  characteristics. 


VP0  V0 


The  circuit  is  modeled  as  in  figure  2a,  where  Cq  is  the  gate 
capacitance  of  T^  and  the  line,  of  length  l,  has  resistance  r  and  capacitance 
c  per  unit  of  length.  Transistor  T^  is  modeled  as  a  (variable)  resistor 
Rq,  to  reflect  the  shape  of  the  v^  =  Vq  characteristic  curve.  In  particu¬ 
lar,  we  approximate  the  latter  as  in  figure  2b  with  two  straight  line 
segments,  meeting  at  the  pinch-off  voltage  VpQ;  the  saturated  regime 
is  modeled  by  a  horizontal  segment  starting  at  (vDg»^)g)  * 

the  so-called  ohmic  regime  is  modeled  by  a  segment  passing  through  the  origin.  We 
shall  now  study  the  general  discharge  regime,  and  later  specialize  it  to 
the  two  regimes  defined  above. 
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2.1  General  solution. 


Let  v(x,t)  and  i(x,t)  denote  the  values  of  the  line  voltage  and 
current  at  abscissa  x  and  time  t,  respectively.  From  Ohm's  law  and  the 
definition  of  capacitance  we  obtain 


5v  di  5v 

ax  =  "ri*  ax  “  "c  at  » 


whence 


(1) 


These  are  instances  of  the  classical  diffusion  equation  (or  heat  equation) , 
which  has  been  assiduously  studied  over  the  past  century.  It  seem  natural 
to  suspect  that  we  are  dealing  with  a  standard  textbook  problem.  However, 
our  boundary  conditions  deserve  special  attention. 

We  assume  that  the  initial  conditions  be  provided  by 

v (x, 0)  =  vQ(x),  x  €  [0 ,£]  (2) 

(or,  alternatively,  i(x,0)  =  i^(x))  where  Vq(x)  is  an  arbitrary  function, 
while  the  boundary  conditions  at  x  =  0  and  x  =  Z  are  supplied  by  the  nature 
of  the  devices,  that  is, 

v  (0,t)  = -RQi(0,t), 

c0  =  iU,t), 

Here  is  a  constant. 

It  is  convenient  to  normalize  time  and  distance,  obtaining  the 
normalized  variables 


t  2  0  , 

(3) 

t  s  o  . 

(4) 
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After  introducing  the  adimensional  parameters  p  =  rji/Rg  and  y  =  c£/Cq,  the 
corresponding  equations  for  voltage  V(|,t)  and  current  I(§,t)  become 


rv  =  Sv 

^=■2  9r 


hf-lf 


(l’b) 


V(?,0)  =  v0(Ag),  ?  €  [0,1] 


(2 '  a)  ,  I(S,0)  -i0(£g),  §€  [0,1]  (2'b) 


f“(0,T)  -pV(0,T)  -  0 


~f(l,T)+Y  ||(1,T)  «  0 


2 

(3'a),  ~(0,T)  -p  |“(0,t)  -  0  (  ) 

ag  5 


(4'a),  ||(1,t)+YI(1,t)  -  0  ( 


The  diffusion  equation  is  normally  solved  by  separation  of  variables. 
Considering  the  current,  we  seek  a  general  solution  of  the  form  I(|,r)  = 
S(§)h(r).  Equation  (l'b)  is  thus  equivalent  to  the  two  equations 


+  y2g  =  0, 

dr 


dh  2, 
d7+^  h  "  0 


for  constant  y.  Any  function  of  the  form  (A  cos y§  +  Bsiny|)e  ^  T  is  a 
solution  of  (l'b);  the  constant  y  is  any  of  the  eigenvalues  of  the  problem, 
i.e.,  any  choice  which  satisfies  the  boundary  value  problem  (3'b),(4'b). 
Specifically,  after  some  obvious  algebra,  from  (3'b)  and  (4 ' b )  we  obtain 
the  characteristic  equation 


tgp.  = 


Y  +  p  y  Y  +  P 


The  infinitely  many  solutions  of  (5)  occur  symmetrically  with  respect  to 
0.  Therefore  we  restrict  ourselves  to  y  >  0.  (A  graphical  display  of  the 
solution  set  is  given  in  figure  3).  The  eigenvalues  [y ^ :  i  =  0,1,...} 
are  indexed  so  that  yg  <  y^  <  ...;  note  that  >  (2i-l)~/2  for  i  £  1.  As 
is  well-known,  to  each  y^  there  corresponds  an  eigenfunction  g  (£)  which 


Y  +  0  M> 


Figure  3.  Illustration  of  the  solutions  of  the  characteristic  equation. 


simultaneously  satisfies 


g  V(?)  +  nfg.d)  =  o, 

gV(0)  -  pg^(0)  =  o, 
g!(l)  +  Yg.(l)  =  o. 


Unfortunately,  relation  (3").  which  is  equivalent  to  pg^(0)  +Pigi(0)=0, 
fails  to  realize  the  classical  Sturm- Liouville  condition  [10],  so  that 
g ^ ( t ) :  i  =  1,2,...]  is  not  a  set  of  orthogonal  functions.  However  by 
defining  the  "inner  product"  of  functions  on  [0,1]  in  the  following  uncon¬ 
ventional  way 


«u,v»  u(|)v(;)d;  +  HWvW 

uo 

It  is  easily  shown  that  the  eigenfunctions  can  be  normalized  so  that 


<(si’gj»  =  5ij’ 

where  S„  is  the  Kronecker  symbol.  Since  (3")  applied  to  the  general 
expression  g(t)  =  Acosu:  +  Bsinuf  yields  A  =  -  (p  /u )  B ,  we  have: 


i: 

C 

b  *  * 


I 
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S±  (?)  =  Gi(sin|ii§ 


-  cosn.5) 

M-i  i 


(8) 


where  is  a  constant.  We  can  now  project  —  in  the  sense  of  our  inner 
product  (6)  —  the  initial  condition  I(§,0)  on  the  set  and  obtain 


l±  =<<KI,0),gi(5)>>  , 


whence  the  general  solution  for  the  current  is 


2 

®  “H  -T 

I(I,T)  -  Z  I,g,(?)e 

i=0 


(9) 


2.2  Analysis  of  the  saturated  regime 

As  mentioned  earlier,  in  the  saturated  regime  starting  at  t  =  0, 
capacitor  Cg  and  the  line  are  at  voltage  Vq,  and  the  transistor  is  modeled 
as  a  current  generator  with  current  value  -1^.  The  circuit  is  modeled  as 
in  figure  4.  Therefore,  boundary  conditions  (3)  must  be  replaced  by  the 


Figure  4.  Model  of  the  saturated  regime, 
new  nonhomogeneous  conditions 


i(0,t)  =  -IQ. 

The  current  !(£,-)  can  thus  be  expressed  as 


I(S.t)  =  ln(5)  +  i,  (?,r) 
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where  I^(?,t)  satisfies  homogeneous  boundary  conditions 

I^O.T)  =  0 
'  3I1 

-g^(l.T)  +  yi^i.t)  =  o 


(3 ) 


with  initial  conditions  I^(§,0)  =  -Iq(§),  while  the  stationary  term  Iq(?) 
satisfies  the  boundary  conditions 

’  V°>  =  _Io 

tt(1)  +  YIo(1)  =  °- 

The  latter,  and  equation  (l'b),  immediately  yield 

V?>  ■  -  l>- 

Turning  now  to  I^(f>r),  note  that  condition  (3,M)  implies  Rq  =  or 
equivalently,  p  =  0.  As  a  consequence  the  boundary  conditions  for  the 
eigenfunctions  {g^(|)l  become 

g.(0)  =  0 

g:<i)  +  vg. (i)  -  o 

which  are  of  the  Sturm-Liouville  type.  Indeed,  from  (8)  the  eigen¬ 
functions  become 


^It  should  be  noted  that  the  eigenvalue  Uq  =  0  does  not  yield  a  valid 
eigenfunction. 


sin  2(i . 
i 


s in  U . i , 
i 


(i  =  l,; 


(i) 
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and  they  form  an  orthonormal  set  in  the  conventional  sense.  The  coefficients 


1^  are  therefore  expressed  as 


(where  (  )  denotes  the  conventional  inner  product)  and 


In  conclusion 


I-l  (?»t)  =  Z  I.g,(?)e 
1  i=l  x 


i(?»t)  *  i Q(^n  s-d  +  ^  i^cs^ ' 1 


2 

“M"  -  T 


(10) 


We  shall  refer  to  the  two  terms  in  the  right  side  of  (10)  as  the  stationary 
and  transient  terms,  respectively. 

The  expression  of  V(1,t)j  the  voltage  at  the  gate  capacitor  end  of 
the  line,  is  obtained  from  I(1,t)  and  the  capacitor  equation,  as 


V(1,T)  = 


v  _££ii  rTm  ^de  -v  Z 

0  c0  JoI(  0  C0(1+Y)T+  C0  i=1 


g.CD 

h— 

.j 


1  -  e 


2 

•^T 


( 


(Note  the  corrective  factor  rc£  due  to  the  normalization  of  time.)  From 
this,  by  integrating  1(5, r)  along  the  line,  we  obtain 

1 

V(0,T)  =  V (1  ,r )  +  r*  J  I (Tl, t )d TJ. 

0 

From  this  expression  for  V(0,t)  we  can  determine  the  time  jp0  at  which 
V(0,t)  *  Vpg,  t.e.  the  time  at  which  the  regime  changes.  Assuming  that 

Vpo/vo  =  0.8,  by  numerical  evaluation  we  have  ascertained  that  for  Y  £  10 

-8 

at  r  =  Tpf)  the  transient  term  of  (10)  is  all  but  negligible  (<  10  .  1(5, TpQ)) 


PO 

Therefore  in  this  range  of  y,  we  may  safely  assume 


I(=,TP0)=  I0(Y^fI  5  " 

as  the  initial  condition  for  the  current  in  the  ohmic  regime. 


t 
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2.3  Analysis  of  the  ohmic  regime 


In  this  case,  the  phenomenon  is  governed  by  (l'b)  (3 ' b )  (4'b)  with 
initial  condition 


;(?.°)  -  i0<7tr  ?  - 1)' 


(2"b) 


This  expression  is  projected  on  the  basis  of  the  eigenfunctions  (8) 
according  to  the  unconventional  rule  (6),  thereby  obtaining 


*i  *  I0(y+T  Hi  -  V 

where  1L  ,  and  (i  =  0,1,...)  are  easily  computable  functions  of  the 
parameters  y, ,  p,  and  y.  It  follows  that 


r«.T)  - 10  = „<7hHi-Ki>Si«>e  1 
1=0 


(  12) 


From  (— a-rjil)  we  readily  obtain  the  expression  of  V(§,t),  as  follows: 


V(S,r)  =V(0,r)-  rl  j  icn,r>d7) 

0 


-E0I(0,T)-til0iIoC7iTH1-K1) 


2  . 

e  J  S±  CT1)<1T1 


that  is , 


-  Vo  1^t7+tHi-^)£i(5).  1 


(  13  ) 


with 


A  r*2  P  1  0  (2) 

fj  (?)  =  g.-(0)+D  g  CDdTl  =  -"“(cosii  +  J  siny,.?) 

1  1  v  x  u.  1U.,  1 

0  i  i 


It  can  be  shown  that  the  { f  (5) }  are  a  set  of  eigenfunctions  of  the 


general  solution  of  (1 1  a) , (3 1  a) ,  (4 ' a) 


3.  Discussion  and  conclusions 

Expressions  (11)  and  (13),  which  respectively  give  the  voltage  V(§,t) 

in  the  saturated  and  ohmic  regimes,  are  the  objective  of  our  analysis. 

In  any  given  technology  the  ratio  y/p  =  cRq/tCq  is  a  constant;  therefore 

only  one  parameter  describes  the  behavior.  Several  discharge  curves  have 

-2-1  2 

been  plotted  in  figure  5,  for  the  values  of  y  *  10  ,10  ,1,  10,  10  ,  1C 

Taking  as  propagation  delay  the  instant  tRR  for  which  V(l,tRR)  =  VxH  — 
we  have  plotted  in  figure  6  the  relation  between  t  and  y.  On  a  purely 


Figure  5.  Discharge  curves  for  y  =  10  ,10  ,...,10 

The  broken  lines  describe  the  discharge  at 
constant  current. 


10-1Vrn-2 


Figure  6.  The  relation  between  t  and  v. 


qualitative  basis,  at  this  point  it  is  interesting  to  observe  the  following 

-  £+ 

facts  (here  we  assume  that  p  £  10  y,  as  in  current  technology): 

(i)  For  small  values  of  y  (roughly,  <  10^),  the  propagation  delay 

is  practically  constant  and  is  determined  by  the  characteristics 
of  the  devices.  This  is  the  unchallenged  domain  of  the  constant 
delay  (or  synchronous  model). 

0  3 

(ii)  For  larger  values  of  y  (roughly,  10  <  y  <  10  ),  the  propagation 

delay  is  basically  proportional  to  y  (i.e.  to  the  wire  length  £). 

This  is  the  domain  of  the  capacitive  model. 

3 

(iii)  For  very  large  values  of  y  (roughly,  y  >  10  ),  the  dependency  of 
the  propagation  delay  upon  y  begins  to  deviate  from  linearity, 
i.e.  the  effects  of  the  dispersive  transmission  line  begin  to  be 
felt.  This  is  the  domain  of  the  diffusion  model. 

On  a  less  qualitative  basis,  we  have  examined  expression  (13)  and 
evaluated  V(§,t)  by  summing  at  first  a  very  large  number  of  terms  of  the 
series  in  the  right  side,  and  next  restricting  the  calculation  to  the  first 
term  (corresponding  to  i  »  0).  Since  Uq  is  very  close  to  0  and  >  (2i-l)rr/2, 
as  was  to  be  expected  the  sum  of  all  other  terms  is  negligible  with  respect 
to  the  first  term.  Therefore  we  shall  now  consider  the  approximate  —  but 
basically  valid  —  expression 

2 

V(S.T)  -  RoVy^T  Ho‘W5)e  °*  (14) 

The  time  constant  of  the  discharge,  in  unnormalized  time  t,  has  the 
expression 


where  —  we  recall  —  Uq  is  the  smallest  positive  solution  of  equation  (5). 
If  in  (3)  we  expand  tgp,Q  in  Taylor  series  we  obtain 
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It  seems  reasonable  to  try  to  define  in  this  diagram  the  regions  of 
validity  of  each  of  the  three  models:  synchronous,  capacitive,  diffusion. 
Specifically,  referring  to  equation  (15),  we  may  (somewhat  arbitrarily) 
define  the  region  of  the  synchronous  model  as  the  one  where  t^  <  2RqCq 
(that  is,  the  time  constant  is  at  most  twice  that  due  to  the  devices  alone); 
by  an  equally  arbitrary  criterion,  we  may  define  the  region  of  the 
capacitive  model  as  the  one  where  e  £  1  (a  deviation  which  at  most  cor¬ 
responds  to  doubling  the  propagation  delay).  This  region  is  shown  unshaded 
in  figure  7.  In  the  same  diagram  each  technology  is  represented  by  a 
straight  line  of  slope  +1,  since  —  as  we  noted  —  in  any  given  technology 
p  =  K^V  (K^ ,  a  constant).  Current  MOS  technology  is  characterized  by  the 
following  parameter  values: 

Feature  width  (\ )  -  2.5  p,m 

Field  oxide  thickness  ■  1  p,m 
Gate  oxide  thickness  =  600  A 
Aluminum  thickness  =  1  |j,m 
Power  supply  voltage  *  5  V 

We  assume  that  2X  and  3 \  be,  respectively,  the  channel  length  of  transistors 

and  the  width  of  aluminum  wires;  in  addition  the  minimum  channel  width  is 

chosen  4\  [1].  Recalling  that  the  resistivity  of  aluminum  is  0.28x10  ^  Qm 

that  the  dielectric  constant  of  is  0.46  X  10~10  F/m,  and  that  the 

-1  2 

electron  mobility  in  Si  is  about  0.8x10  m  /Vsec  (we  refer  here  to  the 
n-channel  portion  of  CMOS),  we  obtain  the  following  values  (see  [1]  [11]): 

IQ  =  0.98  mA,  Rq  =  4.05X103  Cl,  CQ  =  4.12X10-2  pF 

3  -10 

r  =  3.78xl0J  fi/m,  c  =  3.46x10  F/m 


whence  p/y  *  rCg/cRg  *  1.10x10  4.  The  corresponding  straight-line  is  shown 
in  figure  7.  In  addition,  assuming  a  maximum  chip  width  of  10  nm,  we  have 
V  34.  The  corresponding  point  is  also  shown  in  figure  7. 

In  a  scaled-down  technology  of  the  foreseeable  future,  not  all 
parameters  are  likely  to  be  changed  according  to  a  fixed  ratio.  Indeed 
it  appears  that  "feature  size",  gate  oxide  thickness,  and  power  supply  will 
be  scaled  down,  while  there  is  a  strong  interest  in  maintaining  the  thick¬ 
nesses  of  both  aluminum  and  field  oxide.  Therefore  a  reasonable  set  of 
parameters  of  a  future  scaled-down  technology  will  be 

Feature  width  =  0.5  p,m 

Field  oxide  thickness  s  1  nm 

o 

Gate  oxide  thickness  =  150  A 

Aluminum  thickness  =  1  p.m 

Power  supply  voltage  =  3  V 

3  -3 

Correspondingly  we  obtain:  Iq  =  1.4  mA,  Rq  =  1.69 X  10  fi,  Cq  =  6.5 X  10  pF 

r  =  1.89  X  104  T./m,  c  =  6.93X  lO*11  F/m,  whence  p/y  =  0.992  X  10"3.  Moreover 

2 

assuming  a  maximum  chip  width  of  50  mm,  we  obtain  y  s:  5.65 X  10  .  The 

max 

corresponding  curve  and  point  are  also  plotted  in  figure  7. 

The  conclusion  we  extract  from  the  preceding  analysis  is  that  not  only 
the  current  but  also  the  projected  MOS-FET  VLSI  technologies  fall  in  the 
domains  of  either  the  synchronous  or  the  capacitive  models.  In  the  latter 

propagation  delay  is  proportional  to  the  length  of  the  wires.  Note,  how¬ 

ever,  that  this  propagation  delay  is  computed  in  the  hypothesis  that  both 
the  driving  and  the  driven  transistors  be  standard  (i.e.,  of  minimum  size). 
However,  by  raising  the  channel  width  of  the  driving  transistor,  the 
current  Iq  increases  and  t  decreases.  Indeed  —  as  suggested  by  Carver- 
Mead  [  1  ]  and  Thompson  [  3  ]  —  if  the  channel  width  is  proportional  to  the 
capacitive  load  for  all  transistors,  one  approaches  constant  propagation 
time  and,  presumably,  current  density  becomes  the  limiting  factor. 


It  must  be  noted,  however,  that  projected  future  technology  may  reach 
the  (conventional)  boundary  of  the  capacitive  model  region.  Beyond  this 
boundary,  a  possible  design  philosophy  —  as  suggested  by  Chazelle-Monier 
[4]  —  is  to  introduce  repeaters  on  long  wires  in  order  to  achieve  delay 
proportional  to  the  wire  length.  Note,  however,  that  in  this  case  we  can 
no  longer  avail  ourselves  of  channel  width  control.  Another  alternative, 
entirely  in  the  realm  of  speculation,  could  be  the  development  of  integrated 
nondispersive  transmission  lines,  where  speed  of  light  considerations  are 
the  controlling  phenomena. 
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